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Abstract 

We show how the symmetry of attractors of equivariant dynamical systems can be 
observed by equivariant projections of the phase space. Equivariant projections have 
long been used, but they can give misleading results if used improperly and have been 
considered untrustworthy. We hnd conditions under which an equivariant projection 
generically shows the correct symmetry of the attractor. 
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1 Introduction 

Dynamical systems with symmetry often have a large-dimensional phase space, which makes 
it difficult to visualize the attractors. For example, the phase space of 6 coupled van der Pol 
oscillators is This paper describes how to choose a low-dimensional projection of the 
phase space that shows the symmetry of a generic attractor. In the example of 6 coupled 
oscillators, a certain 3-dimensional projection generically shows the correct symmetry of the 
attractor. 

Attractors of an equivariant dynamical system need not have the full symmetry of the 
system. For example, if one tries to balance a ruler on its end, it will inevitably fall either to 
the right or to the left. Even though the system has Z 2 symmetry, each of the two attractors 
(right and left) has trivial symmetry. The initial conditions determine which of the two 
attractors is chosen. This is sometimes called spontaneous symmetry breaking. 

Attractors have two different types of symmetry: instantaneous symmetry which holds at 
each time step, and average symmetry which is a symmetry of the attractor as a whole. It is 
quite easy to identify the instantaneous symmetry from the coordinates of one point on the 
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attractor. Therefore, we concentrate on finding the average symmetry of an attractor when 
its instantaneous symmetry is known. This allows lower-dimensional projections compared 
to an observation space that can distinguish any possible symmetry. Our theory also suggests 
projections to use in physical experiments, where the symmetry of the system is approximate. 

Our main result is that a projection to a representation space generically shows the proper 
symmetry of the attractor if the instantaneous symmetry of the attractor is an isotropy 
subgroup of the representation space. We must say “generically” because it is always possible 
to choose a set whose projection has too much symmetry. We prove that an open dense set of 
equivariant perturbations of the attractor project down to a set with the correct symmetry. 
This is what we mean by “generically.” See Theorem |2.6| for details. 


Quite often, the attractor has trivial instantaneous symmetry, and Corollary p..l6| applies. 
In this case a faithful representation space generically shows the correct symmetry of the 
attractor. It is easy to see from the projection if the attractor does indeed have points 
with trivial isotropy. We recommend that a faithful representation space be used to explore 
a symmetric dynamical system. If an attractor with nontrivial instantaneous symmetry is 
observed, then choose a new projection that is “tailored” to the instantaneous symmetry of 
the attractor. 

For example, in Section 3 we consider 6 oscillators coupled in a ring with Dg x Z 2 
symmetry. The Z 2 symmetry is included because each term in the ordinary differential 
equation (ODE) is odd. A properly chosen equivariant projection of onto shows the 
correct symmetry of a generic attractor, provided the attractor contains points with trivial 
symmetry (that is, the attractor has trivial instantaneous symmetry.) Three-dimensional 
space is still rather difficult to visualize. Fortunately the symmetry can be determined by 
a pair of two-dimensional projections of R^: a “top view” and a “side view,” as shown in 
Figure ||. Roughly speaking, the top view shows the Dg symmetry, and the side view is 
needed to distinguish between the inversion of Z 2 and the 180° rotation in Dg. 

When the attractor in the 6 oscillator system has nontrivial instantaneous symmetry, 
then we can then choose a specially tailored one- or two-dimensional projection that shows 
the symmetry of the attractor. (See Figures and ^) On the other hand, we would need a 
6-dimensional observation space to determine the symmetry of an attractor with any possible 
instantaneous symmetry. 

The recent paper “Detecting the Symmetry of Attractors,” by Barany et al. |Q, intro¬ 
duces a method of determining the symmetry of an attractor by averaging over the attractor 
to obtain a point in high-dimensional representation space. The symmetry of the attractor 
is the same (generically) as the symmetry of the point. 

Several papers |1^, |[ have developed the theory of symmetry detectives, and an al¬ 
ternative approach is described in |]A]. Our paper is inspired by “Detecting the Symmetry 
of Attractors from Observations,” by Ashwin & Nicol [0. They separate the process into 
two steps: First the attractor is mapped, via an equivariant map, to an isotropy equiva¬ 
lent observation space, then another equivariant map to a distinguishing representation is 
integrated over the image of the attractor. This idea was already present in the original 
detectives paper |^, where Barany et al. discuss patterns in a square domain. 

We argue that the observation step is all one needs in many cases. The symmetry of 
the attractor can easily be observed visually in a projection of the attractor. The averaging 
(detective) step converts this beautiful picture to a point, and leads to a table of numbers. 
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If a number is zero, the attractor has a certain symmetry. The number is never exactly 
zero, so cutoffs must be chosen. This is not an easy task, especially since the appropriate 
cutoff depends on the dynamics of the attractor. On the other hand, the human brain is 
excellent at pattern recognition. The correct symmetry can be identihed visually in cases 
where the averaging method is inconclusive or misleading. For example, near a bifurcation 
where two conjugate attractors “glue” together, the eye can recognize when the trajectory 
“jumps” from one side to another. Then we know that the attractor has the 2-fold symmetry, 
although the detective method would not identify that symmetry unless an extremely long 
trajectory were averaged. 

The biggest advantage of the detective method is that the process can be automated 
so that the computer can produce scans of parameter space showing the symmetry of the 
attractor(s). This is an extremely useful tool for exploring a one- or two-parameter family. 
Given the right software and a powerful computer, the equivariant projection technique could 
be used to make a movie of a one-parameter family of attractors. 

The remainder of the paper is organized as follows: In Section 1 we give some background 
in group representation theory and describe our results. In Section 2 we prove the theorems. 
In Section 3, examples of 5 and 6 coupled van der Pol oscillators demonstrate our general 
theory. 


1.1 Group actions and group representations 


Throughout this paper, we let T denote a hnite group. We say that a group T acts on a set 
M if there is a map T x M —>■ M : ( 7 , x) h-^ 7 ■ a; which satishes 7 ■ (a ■ x) = ( 7 ( 7 ) • x for all 
7 and a in T, and 1 ■ x = x, where 1 is the identity element in T. (We use the “•” for the 
group action and concatenation for group multiplication.) If 1/ is a (real) vector space, and 
the action of T onV is linear, we say that G is a (real) representation space for T . (We will 
sometimes say “G is a representation of T,” even though the homomorphism T — GL{V) 
induced by the action is the representation.) See [|T^ ch. 4] or |T^ for details on group 
actions. 


Definition 1.1 A group V acts faithfully onV if ■ x = x for all x E V only i /7 = 1. 

In other words, the group action is faithful if no nontrivial element in the group acts trivially. 

Definition 1.2 Let T act on V and on W. Then a function f : V ^ W is T-equivariant if 
/(^ . x) = 7 • /(x) for all 7 G T. 

That is, the T action commutes with /. Note that the two “-’’s in the dehnition of equiv- 
ariance are (possibly) different group actions. Since we know that / maps V to W, this 
notation is not ambiguous. In fact, we will often drop the “•” altogether. 

Example 1.3 Let the dihedral group D 3 = {p.n \ = (p^)^ = 1) act on as 

follows: 

P- (Xo,Xi,X 2 ) = (x 2 ,Xo,Xi) , K ■ (Xo,Xi,X 2 ) = (Xo,X 2 ,Xi). (1) 

This action consists of permutations of the coordinates; p is a 120° rotation about the diagonal 
Xo = Xi = X 2 and k is a reflection through the plane Xi = X 2 . We have presented the group 
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Dg in terms of generators and relations, and the group action is uniguely defined by the 
action of the generators. An example of an eguivariant function is 

/ : ^ R^; (xq, Xi, X 2 ) ^ {xoxfxl, xix^Xq, X2xlxl), (2) 

where the D 3 actions on the domain and the target space are defined in Eguation &■ Another 
eguivariant function (more specifically an invariant function) is 

: R^ ^ R; (xo,xi,X2) xo + xi + X2, (3) 

where D 3 acts trivially on the target space of g. 


1.2 Irreducible Representations 

A group action on V induces a natural action on subsets U CV, namely 'fU := {y-x \ x E U}. 
We say that U is T -invariant if yt/ = U for all 7 G F. We are most interested in F-invariant 
sets that are also subspaces of V. 

Definition 1.4 A representation space V is F-irreducible if it has no T-invariant subspaces 
other than 0 and V. A representation which is not irreducible is called F-reducible 


Theorem 1.5 (Complete Reducibility) If V is a finite group acting linearly on a finite di¬ 
mensional vector space V, then V decomposes as 

1/ = W © V 2 © • • • © W (4) 


where each of the Vj is F -irreducible. 


A proof is found in |]K 


p. 7]. The complete reducibility theorem says that every represen¬ 
tation is a direct sum of irreducible representations. Recall that V = Ri©V 2 means that every 
element in V can be written uniguely as the sum of a vector in Ifi and a vector in V 2 - The 
direct sum is the Cartesian product of two vector spaces which inherits vector space addition 
and scalar multiplication in “the natural” way, i.e. (xi, X 2 ) 4-(wi, 1 ^ 2 ) = (vi-hwi, V 2 + UJ 2 ) and 
n(xi, X 2 ) = (avi, av 2 ). Thus, the notation R^ = R © R is more explicit than R^ = R x R. 

As a consequence of complete reducibility, we can choose a basis in which the matrices of 
the group representation are block-diagonal, with each block corresponding to an irreducible 
representation space. Furthermore, the basis can be chosen so that the matrices in the group 
representation are orthogonal. 

One can always obtain an equivariant map by linear projection onto an invariant sub¬ 
space. 


Example 1.6 Let D 3 act on R^ as in Example |i. A| . The decomposition into irreducible 
subspaces is R^ = Vq © hi, shown in Eigure [|. The projections onto the two irreducible 
subspaces are 0o : R^ —^ Ro — R : R^ —>• Vi = C = R^, defined by 


0o(x) = Xo + Xi + X 2 . 

and 

The D 3 action on Vq is trivial, and D 3 acts on Vi as follows: 

p ■ z = , K ■ z = z. 


(5) 

( 6 ) 

(7) 
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Figure 1 : 
subspaces: 


The permutation action of D 3 on R^. (a) The projections onto the irreducible 
= Vo © Fi. (b) The hxed point subspaces, described in Section \LA . 


1.3 Dynamical systems 

Let F be a finite-dimensional real vector space. There are two ways to construct a dynamical 
system from a map / : F —> F: 

1. The function can be used to generate a discrete system in which the orbits, 

f f f f /Q\ 

are considered for various starting values xq G F. 

2. The function can be used to generate the vector field of a first-order ordinary differential 
equation (ODE): 

X = f{x) (9) 

In this case, the orbits are the solution curves to the ODE. 

If / dehnes a dynamical system in either of these ways, we will say that the dynamical 
system stems from f. 

For either a discrete system or an ODE, attractors are periodic, quasi-periodic, or chaotic 
orbits to which nearby orbits ‘tend.’ To define an attractor, we first define the cu-limit set 
of a point x E V, denoted uj{x), to be all limit points of the forward orbit through x. There 
is no universally accepted dehnition of an attractor, but we take the following: 

Definition 1.7 A set A is an attractor for a dynamical system if all of the following hold: 

1. A is Liapunov stable. (Roughly, if xq is close to A, then Xn is close to A for all n > 0.) 

2. There is an open neighborhood U of A such that A contains uj{x) for all x eU. 

3. A contains a dense orbit. That is, A = uj{x) for some x E A. 
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As a consequence of this definition, an attractor is closed and dynamically invariant. Fur¬ 
thermore, an attractor has an open basin of attraction. 

In this paper we are concerned with equivariant dynamical systems. In other words, the 
dynamical systems stem from equivariant functions f ■. V ^ V. We assume that the F action 
on V is faithful; otherwise we could factor out the kernel of the group action. 


Lemma 1.8 If A is an attractor for a dynamical system which stems from a T-equivariant 
function, then for each element'j ofT, 

'yA = A or jA P\ A = fli. (10) 


Lemma |1.8| is proved in [g]. The idea is that if A is an attractor, then 7 A must also be 
an attractor. Since there is a dense orbit in A, and since A is contained within an open 
basin of attraction, jA is either identical to A or disjoint from A. This dichotomy helps 


the correlation method of |T^ work quite well, since the correlation measures the extent to 
which “fattened” versions of A and 7 A overlap. 

The open basin property of attractors is crucial in the proof Lemma |^. Other dehnitions 
of attractors do not lead to the dichotomy 7 A = A or 7 A n A = 0. For example, attractors 
with riddled basins [|I| or those stuck onto an invariant subspace P] do not have open basins 
of attraction, and they do not satisfy the dichotomy. Without an open basin of attraction an 
attractor is destroyed by the addition of noise. In physical applications a properly dehned 
attractor should have an open basin of attraction. 


Lemma l.S leads to the dehnition 


:= {A C 1/ I A is compact and 7 A fl A = 0 or 7 A = A V 7 G F}. (11) 

We think of AiV) as the set of all possible attractors of an equivariant map f : V ^ V. 


1.4 Isotropy Subgroups and Fixed Point Subspaces 

Assume that F acts linearly on V. The isotropy of a single point x G F is 

S(x) = {7 G F I 7 x = x}. (12) 

It is easy to see that S(x) is a subgroup of F. We leave it as an exercise to the reader to 
show that F acts faithfully on V if, and only if, V has a point with trivial isotropy. (Recall 
that F is hnite.) 

If H = S(x) for some x G R then H is said to be an isotropy subgroup for the action of 
F on V. In general, not all subgroups of F are isotropy subgroups for a given representation 
space. For example, consider the D 3 action on dehned in Example |^. The subgroup 
Z 3 := (p) is not an isotropy subgroup, because all points which are hxed by Z 3 are on 
the diagonal, and these points have isotropy D 3 . (When we write subgroups in terms of 
generators, the relations are “inherited” from the parent group. Hence, (p) < D 3 is dehned 
to be (p I = 1), since the relation p^ = 1 is in the dehnition of D 3 given in Example |1.3| .) 

Ashwin & Nicol introduce the notion of isotropy equivalence, which is central to our 
work. 
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Definition 1.9 Two representation spaces Vi and V 2 of a finite group T are isotropy equiv¬ 
alent, written as Vi ~ V 2 , if H is an isotropy subgroup for the action ofV on Vi if and only 
if H is an isotropy subgroup for the action ofV on V 2 . 

In other words, Vi ~ V 2 when the two representation spaces have precisely the same isotropy 
subgroups. In section |] we present several lemmas concerning isotropy equivalence. 

Example 1.10 Recall the D 3 action on = Vq © Ih described in Examples 0 and 
Since Vi and have the same isotropy subgroups, namely all subgroups 0 /D 3 except Z 3 := 
(p), they are isotropy eguivalent: Vi rs-/ 

For ^4 C 1/ we define two types of symmetry: 

1 . The instantaneous symmetry of A is T{A) = {7 G F | 7 a: = a: Va: G A]. This is also 
called the pointwise symmetry, or isotropy, of A. 

2. The average symmetry of A is S(y4) = {7 G F | 7^4 = A]. This is also called the 
setwise symmetry of A. 

We will usually think of A as an attractor, and the names “instantaneous” and “average” 
reflect this. For singleton sets, S({a:}) = T({x}) = S(x). 

The instantaneous symmetry of a set is always an isotropy subgroup. To see this, note 
that 

T{A) = n SW- (13) 

x&A 

In Lemma |2.3| we show that the intersection of two isotropy subgroups is an isotropy sub¬ 
group. Thus T{A) is an isotropy subgroup. (Since the group F is finite, we need only consider 
finite intersections.) 

While T{A) must be an isotropy subgroup, Fl(y4) can be any subgroup of F, provided the 
group action is faithful. To see this, let x be a point in V with trivial isotropy, and let JT be 
any subgroup of F. Define Hx := {hx \ h G FT}. Then S(iFx) = H. For example, consider 
the D 3 action on R^ defined in Example The set Z 3 ( 0 ,1, 2) = {(0,1, 2), (1, 2, 0), (2, 0,1)} 
has average symmetry Z 3 . (Recall that Z 3 is not an isotropy subgroup for this action.) 

For F acting linearly on V, we define the fixed point subspace of FT < F to be: 

Fixy(FF) = {x G R I h-x = x Vh G Ff} (14) 

We summarize some facts mentioned previously and list some basic properties of fixed point 
subspaces and isotropy subgroups for a linear action of F on R: 

1. Fixy(Ff) is a vector space for any FT < F. 

2. T{A) is an isotropy subgroup for any R C R. 

3. If FF < G then Fixy(G) C Fixy(FF). 

4. IfACB then T{B) < T{A). 

5. Fixi/(T( 74 )) = A ii and only if R is a fixed point subspace. 


1.6 
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6 . T(Fixv'(-f:r)) = H ii and only if H is an isotropy snbgronp. 

Thns, there is an inclnsion-reversing isomorphism between isotropy snbgronps and hxed 
point snbspaces. 

To simplify the notation, we dehne Fixv'( 7 ) := Fixv'(( 7 )) for a single element 7 G F. 
Also, the snbscript “Id” is often omitted. Fignre |l|(b) shows the hxed point spaces of the 
D 3 action on dehned in Example |1.3| . Since k interchanges xi and X 2 we have Fix(K) = 
{x G I Xi = X 2 ]- The diagonal in R^ is Fix(Z 3 ) = Fix(D 3 ). As an example of property 
6 , note that Z 3 is not an isotropy snbgronp of D 3 , and T(Fix(Z 3 )) = 037 ^ Z 3 . 

Fixed point snbspaces are not necessarily F-invariant, bnt they are very important in 
eqnivariant dynamical systems becanse of the following lemma: 

Lemma 1.11 Let F act linearly on V and let f : V ^ V be T-equivariant. IfW^LV is a 
fixed point subspace, then f{W) C W. 

Proof: Let H be the isotropy snbgronp for which W = Fix(iL). Then, for each x G IF and 
each h e H, hf{x) = f{hx) = /(x). Hence /(x) G IF. □ 


If a dynamical system stems from an eqnivariant fnnction /, Lemma p..ll| implies that 
hxed point spaces are dynamically invariant. If the dynamics is an ODE, and the hxed 
point space has codimension one, then the phase space is divided into dynamically invariant 
sectors. For example, consider a 3-dimensional ODE with D 3 symmetry. Fignre |I](b) shows 
that there are 6 wedge-shaped open regions ont of which a trajectory cannot escape. (For 
discrete dynamics in R^, an orbit can hop from one wedge to another.) 

For three conpled oscillators with D 3 symmetry the phase space is R®, one copy of R^ 
for position and another copy for velocity. Then FixR 6 (K) is codimension-two and there are 
no trapped open regions in phase space. In other words, the set of points in R® with trivial 
isotropy is connected. 

It is now easy to see that points with isotropy T{A) are dense in the attractor A\ Any 
dense orbit must contain only points with isotropy precisely T{A), otherwise it would be 
“stuck” in a different hxed point subspace. However, there can be points in A with larger 
isotropy if they are accumulation points of a dense orbit. 

On the other hand, there are sets U where no points in U have isotropy T{U). For 
example, let U be the union of the three hxed point subspaces shown in Figure |l|(b). No 
points in U have trivial isotropy, but T{U) = 1. 


1.5 Invariant Poincare Sections 

Note that F-invariant subspaces of V are not necessarily invariant under the dynamics. 
However, F-invariant subsets are important since they are candidates for invariant Poincare 
sections, as stated in the next theorem. Consider, for example, an eqnivariant ODE with the 
group action shown in Figure |^. The F-invariant subspace Fi is not necessarily dynamically 
invariant, but it is an invariant Poincare section: The symmetry of the attractor A C R^ is 
the same as the symmetry of the intersection of A and F, provided the intersection is not 
empty. 





Theorem 1.12 Ashwin & Nicol Lemma 6.2], Suppose A is an attractor for f : V ^ V, 
and P <ZV satisfies (a) ■yP = P for all y eT and (b) a dense orbit in A intersects P. (Such 
a P is called an invariant Poincare section.j Then T{A (1 P) = 2(^4). 


The proof, fonnd in [Q], is qnite easy considering that this important theorem seems to be 
new. 

Note that P can be any P-invariant set. It can have any dimension, and it need not be 
a linear snbspace. Furthermore, P is not really a Poincare section. There is no requirement 
that the flow is transverse to P, and a point is plotted every time the trajectory hits P, not 
just when it crosses in one direction. 

We will only consider invariant Poincare sections P which are P-invariant codimension- 
one subspaces. These are easily classified: they are the orthogonal complements to the 


one-dimensional representation spaces in the decomposition of V described in Theorem ^ 
For example, with D 3 acting on as in Example 


and Figure |l|, the unique codimension- 
one invariant subspace is Xf, which is the solution to fio^x) = Xq + Xi + X 2 = X). For any 
linear group action on V, if is the linear projection onto a one-dimensional irreducible 
representation of V, then (j){x) = 0 gives the equation for a codimension-one F-invariant 
subspace. 


1.6 Projections which Generically Preserve Symmetry 

In this section we present what is needed to apply our theory to the observation of symmetries 
of attractors. Those who are interested can read the theorems and proofs in Section 2. Others 
can skip Section 2 and go directly to section 3 where these ideas are applied to systems of 
coupled oscillators. 

A crucial observation is that an equivariant map can never decrease the symmetry of a 
set: 

Lemma 1.13 Suppose that F acts linearly on M and V. Then for any equivariant map 
cj) : M ^ V, and any A C M, 

S(A) < S(0(A)) and T{A) < T{(j){A)). (15) 

Proof: Let p G E(A). Then p4>{A) = 4>{pA) = (j){A). Thus p G S((/)(A)). The proof for 
T{A) is similar. □ 


Note that if the map (j) is invertible we have equality in Lemma |1.13| , but in general we 
are interested in non-invertible maps for which the domain is a high-dimensional phase space 
and the range is a low-dimensional observation space. 

The following example shows two ways in which an equivariant map can increase the 
symmetry of a set. 


Example 1.14 Let Z 2 = (k | = 1) act on by k ■ {x,y) = {—x,y). Consider the two 
equivariant projections fifix^y) = x and (f) 2 {x^y) = y. For any set A C R^, S( 02 ( 2 l)) = 
r(02(2l)) = Z 2 since the F action on the y-axis is trivial. The projection onto the faithful 
representation can also lead to the wrong symmetry. For example let A = {(—1,0), (1,1)}. 
Then S(A) = 1, but E(0i(A)) = S({-1,1}) = Z 2 . 


9 




This example shows that we cannot observe the correct symmetry of all sets. The best 
we can hope for is that S(0(y4)) = S(y4) when A is a “generic” set within its symmetry class. 
To make this precise, we imagine pertnrbing the set A to i^iA), where V’ is an eqnivariant 
map from M to M. (If we did not reqnire V’ to be eqnivariant, then “most” pertnrbations 
wonld give a set with trivial symmetry.) The following example shows how an eqnivariant 
pertnrbation of A can get rid of spnrions symmetries in the projection. 


Example 1.15 Consider the Z 2 action on defined in Example 1.14 - Define the equiv- 
ariant map y) = {x + ex sin ny, y). For the set A defined in Example the perturbed 
set is 'fi{A) = {(—1, 0), (1 + e, 1)}, so fiiifij^A)) = { — 1,1 + e}, which does have the same 
(trivial) symmetry as A for e 7 ^ 0. 


To make progress, we make the additional assnmption that A satisfies the dichotomy de¬ 
scribed in Lemma |1.8| . Figure |] shows what can happen withont this assnmption. Let Z 2 
act on as described in Example |1.14 The set A is made np of two arcs. The set B, 
consisting of the 4 endpoints of the arcs, has E{B) = Z 2 , even thongh S(y4) = 1. (Hence 
A does not satisfy the dichotomy.) The projection fifiA) is two intervals, the endpoints of 
which are 0i(R), so S(0i(H)) = Z 2 7 ^ = 1- For any eqnivariant ip : R^ —> R^ we 

have S(0i oip[B)) = Z 2 , since eqnivariant maps cannot decrease the symmetry. If ip is near 
the identity map, then pi o ip{A) is two line segments with endpoints pi o p{B). Hence 
S(0i o p{A)) = Z 2 , for an open set of p, even thongh S(H) = 1. 





B 


t t . 


Fignre 2 : The projection of the set A onto the a:-axis has S = Z 2 , even thongh A has 
trivial symmetry. Any near-identity eqnivariant pertnrbation of A will project down to a 
symmetric set becanse the endpoints, denoted by R, have Z 2 symmetry. 


It is fortnnate that attractors satisfy the dichotomy of Lemma |^, since in this case the 
problem shown in Fignre ^ does not occnr. Recall that A{M) is the set of all compact snbsets 
of M which satisfy the dichotomy. Ashwin & Nicol 0] prove that E(A) = S(<;/>(A)) for all p in 
an open dense set of eqnivariant maps from M to V, provided that V is isotropy eqnivalent 
to M, A G A{M), and points in A with isotropy T{A) are dense in A. However, this theorem 
does not imply that an eqnivariant linear map p will preserve the symmetry, and these are 
the ones we feel are most natnral for observing the attractor. We take a different approach, 
and pertnrb A by an eqnivariant map p : M ^ M as in Example |1.15| . Let Cf{M) be the set 
of /c-times continnonsly differentiable F-eqnivariant fnnctions from M to M. Onr theorems 
give conditions on linear maps p : M ^ V which ensnre that S(0 o p{A)) = S(A) for all 
p & O, where O is open and dense in Cp(M). 

We now state a corollary of our main theorem, which is stated and proved in section 2. 
Consider the case where the attractor has trivial isotropy, that is T{A) = 1. Note that if any 
points in the projection p{A) have trivial isotropy, then T(A) = 1 since T{A) < T{p{A)) = 1. 
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Corollary 1.16 Let M and V be finite dimensional real representation spaces of a finite 
group r. Let (p be a T-equivariant linear map from M onto V. Suppose that A G A{M), 
and assume that V is a faithful representation ofT. If there is a point a E A such that 
T{(j){a)) = 1, then 

T{A) = 1 and S(/l) = E(0 o 'fi{A)) (16) 

for all E O, where O is an open dense subset of Cf{M). 

In practice, this means that the symmetry of the projected attractor 4>{A) is generically 
the same as the symmetry of A. By “generically,” we mean that if the symmetry of the 
projection is larger than the true symmetry of A, then “most” equivariant perturbations 
V’(^) will project down to a set with the correct symmetry. Since attractors of dynamical 
systems are quite “generic” sets, we can be conhdent that 0 shows the true symmetry if V 
is a faithful representation. 

To observe the symmetry of the attractor with nontrivial instantaneous symmetry, we 
need an observation space V for which T{A) is an isotropy group, as described in The¬ 
orem \2.6[ This V is often smaller than a faithful representation space. Attractors with 
T{A) 7 ^ 1 are observed in the example of 6 coupled van der Pol oscillators in Section |3.3| . 

2 Theory 

2.1 Isotropy Equivalence 

In this section we develop some theory regarding the notion of isotropy equivalent repre¬ 
sentations, dehned in Equation (|L^). We use the notation Vi ~ V 2 to indicate isotropy 
equivalence when the group T is understood. Note that ~ is an equivalence relation. 

The following three lemmas develop several properties of isotropy subgroups. 

Lemma 2.1 Let T be a finite group with representation V. Let H < T. The following are 
equivalent: 

1. H is an isotropy subgroup. 

2. For all G <T such that H <G, Fixv'(G) = Fixy(il) if and only if G=H. 

3. Fixy(G) C Fix\/(if) for all G <T such that H < G. 

We take C and < to mean proper subset and proper subgroup respectively. 

Proof: 

(1 ^ 2) Assume that H is an isotropy subgroup. Let G < T such that H < G. Clearly 
a G = H then Fixy(G) = Fixv{H). Suppose Fixv'(G) = Fixy(iL). Since H is an 
isotropy subgroup, there exists v E Fixy(iL) such that H = F(r>). Thus v E Fixy(G) 
and since H = T{v), G < H . Thus G = H. 

(2^3) This is obvious. 
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(3 ^ 1) Assume that Fixy(G) C Fixv{H) for all G < F such that H < G. Let Q = 
{G < F I Ff < G}. Then Fixy(G) C Fixy(iF) for all G & Q. Since Q is hnite, 
IJGeeFixy(G) 7 ^ Fixy(iF). Consequently there exists 

V e Fixy(iF) \ y Fixy(G). (17) 

G&g 

Thus Ff = S(n), so iF is an isotropy subgroup. 


□ 


Lemma 2.2 Let T be a finite group with representation V. Given FF < F there exists a 
unique F <r such that H < F, Fixy(FF) = Fixy(F^), and F is an isotropy subgroup. 

Proof: Let F be the subgroup of pointwise symmetries of Fixy(FF), i.e.: 

F = { 7 GF| 7 u = u VuG Fixy(FF)}. (18) 

Clearly, H < F and Fixy(F^) = Fixy(FF). That F is an isotropy subgroup follows from 
Lemma M The uniqueness of F is obvious. □ 


Lemma 2.3 Let FF, G < F be isotropy subgroups for a representation V. Then FF fl G is an 
isotropy subgroup. 

Proof: Let F < F such that G fl FF < F. Then there exists f ^ F such that / ^ G or 
f ^ H. Without loss of generality assume f ^ H. FF is an isotropy subgroup, so there exists 
V eV such that FF = S(u). Thus fv 7^ v. Since G fl FF < FF, u G Fixy(G fl FF). Therefore, 
Fixy(F) C Fixy(G fl FF). By Lemma ^T|, FF fl G is an isotropy subgroup. □ 


Note that F is always an isotropy subgroup. If the action of F on P is faithful, then 
1 is an isotropy subgroup, otherwise the kernel of the action is an isotropy subgroup. In 
fact. Lemma |2.2| indicates that every subgroup is contained in a unique “smallest” isotropy 
subgroup. 

Now we present two lemmas. The hrst of these provides a method for showing that 
two representations are isotropy equivalent. The second allows us to reduce the size of 
representations by removing isotropically equivalent representations which are “redundant.” 


Lemma 2.4 Let T be a finite group with representations Vi and V 2 . Suppose that, given 
H <G<T: 

Fixy, (G) = Fixy, (FF) ^ Fixy, (G) = Fixy, (FF). (19) 


Then Vi ~ V 2 - 


Proof: This is a straight forward application of Lemma 


□ 


Lemma 2.5 Let T be a finite group with representations V, Vi, and V2. F/ Pi ~ V2 then 
p'ePi ©p ~ P'ePi. 
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Proof: Clearly if H is an isotropy subgroup for V © Vi then H is an isotropy subgroup for 
r © Pi © P2. 

Now let H be an isotropy subgroup for P' © Pi © P2. Then H = S((u', ui, V2)) for some 
iy',Vi,V2) G P' © Pi © 1 / 2 - Since the elements of H must £x {y',Vi,V2) coordinate-wise, 
H < F n Gi n (^2 where F = S(n'), Gi = S(ni), G2 = S(r;2)- However, any element of 
F n Gi n G2 hxes (n', ^1,^2), so F fl Gi fl G2 = TT. 

By assumption, G2 is an isotropy subgroup for Pi. Thus, by Lemma Gi fl G2 
is an isotropy subgroup of Pi. Suppose that Gi fl G 2 = S(t') for some v G Pi. Then 
F = F n Gi n G 2 = H(u') n T,{v) = E((u',r;)). Therefore H is an isotropy subgroup for 
P'©Pi. □ 


2.2 Cy Functions 

In what follows, we use the notation ||a;|| for the Euclidean norm of a point x G Given 
a point a E X, where X is a metric space, let B{a, 5) denote the ball around a of radius 6. 
We indicate the closure of a subset F of a topological space by B. 

Given hnite dimensional real vector spaces M and P, let G^{M, V) denote the set of all 
fc-times continuously differentiable functions from M to P, and define G^{M) := G^{M, M). 
Given 0 G G^{M, P), define the standard norm of (j) by 

||0||c-fc := sup{||F*^-^V(^)ll I 0 < j < fc and x G M}. (20) 


When j > 1, DG'>(1)[x') is the derivative of cj) at x. When j = 0, DG'^(f){x) = (j){x). 
For j > 1, D^^'^(p(x) can be defined as a j linear map from M into P. In Equation (|^) , 
||F^-^^0(a;)|| denotes the operator norm of D^^^(p(x). 

If M and P are representations of a finite group T let Gy{M,V) denote the set of G^ 
functions which are T-equivariant. Define G^{M) = G^iM, M). We use the strong topology 


on G^(M, P) and Gf (M, P) fn 


2.3 The Main Theorem 


We feel that equivariant linear functions are the most natural way to visualize the phase 
space. However, it is certainly not true that linear functions preserve the symmetry of every 
subset of M. Our main theorem says that if A satisfies the dichotomy 'jA = H or yHflH = 0, 
and if the linear map (j) : M ^ V satisfies certain conditions, then almost all equivariant 
perturbations of A have the correct symmetry when observed in P. 

Before we prove Theorem p.6| , recall that given any 0 G Gp(M, P) and A C M, E(H) < 


S(0(H)) and T{A) < T{(j){A)) (Lemma |1.13|) . Recall, also, that we defined: 

A{M) = {H C M I H is compact and 7H = H or 7H nH = 0 V7Gr}. 


( 21 ) 


Theorem 2.6 Let M and V be finite dimensional real representation spaces of a finite group 
T. Let (f) be a T-equivariant linear map from M onto V. Suppose that A G A{M) satisfies 

1. T{A) is an isotropy subgroup for V. 
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2. There exists a E A such that S(a) = T{A). 

Then for each positive integer k there exists an open dense set O C C^^M) such that 
S(0(^(a))) = S(a), T(0(^(A))) = T{A), and S(0(^(A))) = S(A) for all fjeO. 

Proof: To simplify notation let C = Cp(M). Define 


0 = {%IjeC\ S(0(^/>(a))) = S(a) and S(0(^(A))) = S(yl)}. (22) 

Then for any E O, T{A) < T{(j){ij{A))) < S(0(^(a))) = S(a) = T{A), so T(0(^(A))) = 
T{A). It remains to show that O is open and dense. 

For each r G T \ S(a), let Wr = {i/j E C \ t E ^(^('^(a)))}. For each p G F \ ^(A), let 
Vp = {fj E C \ p E S(0(^/>(yl)))}. Clearly, 


O 


C\ 


U WrU U Vp), 

TGr\S(a) pGr\S(A) / 


(23) 


We show that each of the sets and Vp are closed and nowhere dense. It follows, since F 
is hnite, that O is open and dense. 

If S(a) = F, the result is trivial. Thus suppose that S(a) < F. 

Let r G F \ S(a). Dehne F : C —R by F{^jJ) = ||r ■ (j){'ip{a)) — (j){'ip{a))\\. Clearly F is 
continuous and Wr = F“^(0), so Wr is closed. To show that the closed set is nowhere 
dense, it suffices to show that the complement of Wt in C is dense. Let e > 0 and fj E W^. 
We must show that there exists E C \ Wr with Wfj — fjWck < e. Since T{A) is an isotropy 
subgroup for V, there exists v E V such that S(n) = S(a) = T{A). Since 0 is onto, there 
exists w E M such that 0(tc) = v. Let 5 > 0 such that B{a,S) fl 'yB{a,6) = 0 for every 
7 G F \ S(a). Such a S exists because F is hnite and A E A. Choose p E C^(M, R) with 
support B{a,S) such that 77 (a) = 1. Note that \\p\\ck < 00 since 77 has compact support. 
Dehne ip E C hj 

iP{x):=ijj{x)+ -w, (24) 

C\\v\\c^ 

where C = I].ygr lib ' 'W^ll- Then 

(l){iP{a)) = (j){i/j{a)) + ^ |S(a)|n. (25) 

<b||hllc"= 


So 

r ■ (piipia)) - (piipia)) = ^ |S(a)|(r ■ n - n) ^ 0. (26) 

Cllhllcfc 

Thus r ^ S(0(7/’(a))), so ip ^ Wr. But 


11 ^ -^llc^ = 


< 


e 

Cllhllc*^ 

e 


^77(7 ^ x)j 

7 Gr 



WvWc^Yl llb-w^ll =e. 

7 Gr 


(27) 
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Now we show that Vp is closed and nowhere dense. If S(yl) = F, we have nothing to 
prove. Therefore assume that S(y4) < F. 

Let p G F\E(y4). Let D{B,C) denote the symmetric Hausdorff distance between compact 
subsets B and C of V. Let F : C — R be dehned by F{7p) = D{p ■ (l){'ijj{A)),ip{A)). Then 
F is continuous and Vp = F“^(0) so Vp is closed. To hnish the proof, we must show that 
R(V’, e) n (C \ Vp) is nonempty for every 'ip E Vp and e > 0. Let p) EVp and e > 0. We have 
already shown that C \ UrWr is open and dense, in particular there exists ip E C such that: 

1- 11^ < e/2. 

2. E(0(^(a))) = S(a). 

If p ^ S(0('^(yl))) then 'ip E B{'ip,e) fl (C \ Vp) and we are done. Thus assume that p E 
Tj{(p{'ip{A))). In other words, (p{'ip{A)) = p ■ (p{'ip{A)) even though An pA = %. We shall 
perturb ip so that an extreme point of the image of y4 in V gets pushed outward from 0 , 
yielding the correct symmetry of the image. See Figure |^. 



Figure 3: The proof in the case where Z 2 acts on by {x,y) h->• {—x,y), and (p{x,y) = x. 
The projection of '0(A) has Z 2 symmetry, even though A does not. We perturb 0 to "0 so that 
0(0(A)) does have the correct symmetry. Note that 0(A) need not satisfy the dichotomy 
that A does. 

Since 0(0(A)) is compact, we can choose to maximal such that t^cppip^a)) E (ppip{A)). 
There exists a' E A such that 0(0(a')) = to(p{%p{a)). Dehne b = 0(0(a')). Note that 6 is a 
boundary point of the image of A in V, since tb ^ (p{'ip{A) for f > 1. Since F acts linearly 
and p E S(0 o0(A)), we conclude that tpb ^ (p{'ip{A)) for f > 1. Further, S(a) = Fl(a') = 
S( 6 ) = T(A), since Fl(a') < E(0(-0(a'))) = E(0('0(a))) = S(a) = T(A) < E(a'). Keep in 
mind that & 0 0, since we are assuming that E(a) < F. 

Let S' = (0 o 0 )“^({tp ■ 6 I 0 < f < 1}) n A. That is, S contains all the points in A whose 
image under 0 o 0 lies in the line segment from 0 to p ■ 6. Since p ^ F1(A), and A G A(M), 
we know that p ■ a' ^ A. Furthermore, 7 a' ^ S for every 7 G F. Since F is hnite, we can 
choose (5 > 0 such that jB^a', <5) fl S' = 0 for every 7 G F and aB{a', 6) fl B{a', 5) = 0 for 
every a ^ S(a'). Choose a non-negative function p G C^(M, R) with support B{a',6) and 
p(a') = 1. Let X denote the characteristic function of the group orbit {■jB{a',6) | 7 G F}. 
Dehne 0 G C by 

H^)-= (l + .,^11^111 II II 7 II (28) 

\ 2 + 0 F|||p||cfe|| 0 x||cfe ^gr / 
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Note that 0('0(x)) = 0(V’(a;)) for every x G S'. Also, (p{'ijj{x)) oc 0('0(x)) for every x G A, 
with a constant of proportionality greater or eqnal to 1. In particular, 


(j){ip{a')) = 1 + 


2 ^+^|r|||hllc^llV’Xllc^ 


|S(aO| = tb. 


(29) 


with f > 1. Given x G A, 0('^(x)) = tp ■ b only if 0 < f < 1 in which case x G S'. Thus, 
p- (j){'ip{a')) ^ 0('^(A)). Therefore V’ ^ Vp. Furthermore, 


||^-^||cfc< 11 ^ 

< 


Ic-fc + 11-0 — V’llc'fc 


< 


... II E ||h(7a^)V’(a:)x(a:)llc'= + e/2 

2^|r| llhllc'=l|t/’Xllc'= + e/2 = e. 


(30) 


2'^+^|r||hilcfc||^/’Xllcfc 


In the second line we used the fact that p{'yx)'ilj{x) = r7(7x)'?/(x)x(x) for all 7 G T. This 
complication is necessary because ||'?/||c'' might be inhnite. In the third line we used the fact, 
easily demonstrated with Leibnitz’ rule, that H/fi'Hcfc < 2 ^||/||(7fc||5f||cfe for any / G G^(M, R) 
and g G 

Thus Vp is closed and nowhere dense. □ 


Note that the two-step procedure of perturbing 'ip io 'ip and then to pj is necessary when 
■^(x) = 0 for all X G M. 


Corollary 1.16 follows from Theorem ^ if we restrict to T(A) = 1, and recall that a 
representation space of a hnite group is faithful if it has points with trivial isotropy. 


3 Oscillators coupled in a ring 


As an example of an equivariant dynamical system, we consider n > 3 coupled van der Pol 
oscillators: 


= Vj - ^(vj-i + yj+i) 

Vj = + %•(« - + PiVj-i + %+i) + 7(%-2 + yj+2) + S{xj_i + Xj+i) 


(31) 


where j G {0,1,..., n — 1} is dehned modulo n. We can think of the oscillators coupled in a 
ring with coupling to their nearest neighbors (the [3 and 5 terms) and their next-to-nearest 
neighbors (the 7 term). We investigated Equation (^I]) using dstool driven by a fourth order 
Runge-Kutta algorithm with hxed time step h = 0.01. 

The equations have the symmetry T = x Z 2 . The symmetry of the regular n-gon, 
D„, is due to the geometric symmetry of the ring. We write the n-fold rotation in as p, 
and the reflection as k. Thus := (p, k : p^ = = (p^)^ = 1). (Note that [lA] and |n 

use C for the rotation.) The symmetry group acts on = R” © R” as two copies of the 
same representation, one for the x coordinates and one for the y coordinates. For simplicity 
we show the action on the x coordinates only: 


p ■ (xo,Xi, . . . ,Xn-l) = (xi,X 2 , . . . ,Xo) 

^ ’ (^0) ^li • • • 1 1 ) (^0) ^n—li ^n—2i ■ ■ ■ i ^l) 


(32) 

(33) 
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Note that k acts trivially when n = 2. This is why we assume n > 3 throughout this paper. 

The Z 2 symmetry in T = x Z 2 is present because every term in Equation (|3T|) has 
odd degree. The equations are symmetric under the inversion through the origin in R^"'; 

■ {x,y) =-ix,y). (34) 

Clearly a commutes with all elements in D„. It is very common for oscillator models to have 
this inversion symmetry. 


3.1 Irreducible Representations of and D„ x Z 2 

When we refer to as the symmetries of a regular n-gon, we are implicitly thinking of 
the action on the plane, which is most simply written in terms of complex coordinates 
zeC^R^: 

p ■ z = , K ■ z = z (35) 

The complex coordinate is used to denote the 2-dimensional real representation T ^0(2); 


/ cos( 27 r/n) — sin( 27 r/n) \ / 1 0 \ 

y sin( 27 r/? 7 .) cos( 27 r/n) y ’ ^ y 0 —1 y 


Equation (^5]) does not give a 1-dimensional comp/ea; representation T —>• [/(I) because of 
the complex conjugate z. 

The group action in Equation (|35|) can be generalized to give an irreducible representation 
space Vm for every integer m, where D„ acts on Zm ^ Vm as follows: 


P ‘ ^ra 



^ ‘ ^m. 


(37) 


The Vm are also irreducible representation spaces of D„ x Z 2 , where the generator of Z 2 acts 
as 

a- Zm = -Zm (38) 

Note that Equation (^5]) is the group action on Vi. The representation space Vm is the same 
as Vn+m and it is conjugate to Vn-m- Hence, we only need consider m in the set 0 < m < n/2. 
When 0 < m < n/2, Vm = C is a (real) 2-dimensional representation space, whereas Vq = R 
and Vn /2 = R (if n is even) are one-dimensional representation spaces. The rotation p G D„ 
acts as a rotation by 27im/n radians in Vm- 

Equation ( |37D does not describe all of the irreducible representations of D„. There is a 
one-dimensional irreducible representations in which p 1 —^ 1 and k 1 —— 1. If n is even, there 
is another one-dimensional irreducible representations in which p 1 —>■ — 1 and k 1 —— 1. 

We now introduce linear equivariant functions of the phase space: Let {x,y) E R^*^, and 
for each integer m dehne Vm ^Vm and Wm E Wm as follows: 

n—1 . n—1 

Vm = 4>m{x,y) := Xj ^ = 'ipmix,y) := Y Vj . (39) 

j=0 j=0 


Figure | gives a graphical interpretation of this sum: The coordinates Xj or yj are mul¬ 
tiplied by various roots of unity in the complex plane. 
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Figure 4: This shows the various roots of unity used to dehne Vm = (f)rn{x,y) and Wm = 
y) for n = 5 and n = 6 oscillators. 


The functions (pm and ipm dehned in Equation (|3^ are equivariant linear projections 
onto the irreducible representation spaces Vm and Wm, respectively, of D„. (These are also 
irreducible representation spaces of D„ x Z 2 .) Note that pm generalizes the projections 0o 
and 01 dehned in Example for n = 3. For a given m, Wm is simply a copy of Vm, and we 
could have written Vm,Wm G Vm- Instead we write Vm G Vm, and Wm ^ Wm to distinguish 
the two copies of the same representation space. Note that pm = Pm+n and pn-m = Pm- We 
obtain a linear change of basis by considering only 0 < m < [n/2j (the hoor of n/2). In 
fact, the decomposition of the phase space described in Lemma is 

Fq © Wi © • • • © W\nl2\ (40) 


M = 

R" © R" = 

Wo © El © • 

■ ■ © V\n/2\ ^ 

action 

is block-diagonalized in 

the (n, w) 1 


pVm — 

Vm 

^ 1 P'^m 







1 



^gi27r/ny 


W. 


-^m 

-w^ 


(41) 


Note that only the irreducible representations of D„ described in Equation ( 0 ) are present 
in the decomposition of the phase space. 

Recall that a linear projection onto a faithful representation space generically preserves 
symmetries of attractors with trivial isotropy (see Corollary p..l 6 | ). 

Lemma 3.1 The projection pi : R^"' Vi is a linear projection onto a faithful representa¬ 
tion space o/D„. For the group D„ x Z 2 , Vi is a faithful representation space if n is odd, 
whereas Vq © Vj is the smallest faithful representation space if n is even. 

Proof: This is straightforward. When the group is D„ x Z 2 and n is even, p'^^^avi = vi, so 
the action on Vi is not faithful. □ 


In order to observe any possible symmetry, including cases with nontrivial instantaneous 
symmetry, we need a representation space that is isotropy equivalent to the phase space 
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Lemma 3.2 The representations Vq and K 0 /D„ or x Z 2 are isomorphic, and hence 
isotropy equivalent, if and only if gcd{q,n) = gcd(r, n). (Recall that gcd(g, n) is the greatest 
common divisor of q and n.) 

Proof: The root of unity in Equation (^) can be replaced by any other primitive root 
of unity, namely where the gcd(p, n) = 1. (Note that p and n are relatively prime 

precisely when gcd(p, n) = 1.) We call the projection cj) when using and 0 when using 

gl27r/n_ 

~ n—1 ^. 

4>m{x, y) = Y1 = (fpmix, y) (42) 

i=o 

But (fm and (fm are isomorphic since they are induced by the outer automorphism of D„ 
dehned hj p ^ p^, k ^ k. Therefore Vm — Vpm provided p is relatively prime to n. (We 
use “=” to denote isomorphism, and to denote isotropy equivalence. Note, however 
that Vm and Vpm are not conjugate representations, and they have different characters.) 
The isomorphic representations correspond to the cycles of the permutation m 1 —>■ pm mod¬ 
ulo n, where gcd(p, n) = 1. The gcd(m, n) is constant for every element of a cycle, and 
different for every cycle. (For example, the cycle decomposition of m 1 —>■ 3m modulo 10 
is (0)(1397)(2684)(5), and the greatest common divisors are 0, 1, 2, and 5, respectively.) 
Hence, Vq = Vj. if gcd(g, n) = gcd(r, n). Conversely, if gcd(g, n) 7 ^ gcd(r, n) the representa¬ 
tions cannot be isotropy equivalent because the general point in Vq has isotropy 
whereas the general point in W has isotropy The proof is unchanged if the group 

is D„ X Z 2 rather than D„. □ 


For example, for n = 10, Vi = V 3 and V 2 = V 4 . 

Theorem 3.3 For the permutation action o/D„ on the smallest representation space 
isotropy equivalent to is the direct sum of all Vm where m divides n and m < n: 

R'" ~ 0 Vjn. (43) 

m|ri., m<n 


ForY)n X Z 2 , the smallest representation space isotropy equivalent to R^"' also includes Vq: 


r2^ ~ 


W © 0 Vm. 

m\n, m<n 


(44) 


Proof: Lemma states that that V ® V ® W ~F'©FifF ~ W . Note that Vm 
and Wm (with the same m) are conjugate and therefore isotropy equivalent, since there is 
an invertible linear change of coordinates which interchanges Xj and yj for all j. Thus, we 
can “throw away” the Wm representations in the decomposition Equation (|40|). Lemma 
implies that we only need Vm where m is a factor of n. (This includes Vn, which is the same 
as Vq-) For Ihe group D„, Vq is the trivial representation and can be eliminated. Finally, we 
cannot get rid of Vq for D„ x Z 2 , since the point ( 1 , 0 ) G Vq © Vi has isotropy D„, but no 
point in Vi has this isotropy (0 G Vq has isotropy D„ x Z 2 ). □ 


For example, for n = 10 coupled oscillators in Equation ( 0 ), the smallest observation 
space which is isotropy equivalent to the phase space R^° is Vq © Vi © V 2 © V 5 , which is 6 
dimensional (l + 2 + 2 + l = 6 ). 
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3.2 Five coupled van der Pol oscillators 

The case of n = 5 van der Pol oscillators is typical of systems with D„ x Z 2 symmetry when 
n is prime. The representations Vi, V 2 , Wi and W 2 are all isotropy eqnivalent and faithfnl. 
(We remark that Vi and V 2 are not lattice eqnivalent, contrary to the claim in |^.) While 
Vo © Vi is the smallest representation space that is isotropy eqnivalent to the 10 -dimensional 
phase space, we only need to look at Vq when the attractor projects to a point at the origin 
in Vi- 

Here is an algorithm to determine the symmetry of the attractor, assnming that H is a 
“generic” set: 

1. If (t)i{A) = 0, then T{A) > D 5 . In words, if the Vi projection of the attractor is a point 
at the origin, then either the oscillators are in-phase (T(A) = D 5 ), or the attractor is 
a point at the origin in ('r(A) = D 5 x Z 2 ). Consider then the projection onto Vq: 

(a) If MA) = 0, then T(A) = E(A) = D 5 x Z 2 . 

(b) If MA) = -MA) + 0, then T{A) = D 5 and S(H) = D 5 x Z 2 . 

(c) If (t>o{A) + -(/.o(d), then T(H) = S(H) = D 5 . 

In practice, we wonld look at the (no, Wq) plane rather than the one dimensional projec¬ 
tion onto Vq. In case (b) the hgnre wonld have 2-foId rotational symmetry in Vq © ITo, 
whereas in case (c) it wonld not. 

2. If 01 (H) 7 ^ 0, then the symmetry of A is the symmetry of the projection onto Vi- 
The instantaneons symmetry is T{A) = Di if the attractor lies in a line of reflection, 
and T{A) = 1 otherwise. Note that D 5 x Z 2 = Dio, and a G S(yl) precisely when 
01 (H) = —0i(H). (Recall that a is the generator of Z 2 .) 

Figure ^ shows the attracting periodic orbit of Equation (|^) for n = 5, a = 2, 0 = —0.5, 
and 7 = 5 = 0. The negative 0 coupling pushes adjacent oscillators out of phase. We show 
the projections onto Vi and the other faithful representations. The large dots in Figure ^ 
are an invariant Poincare section obtained by plotting a point whenever wq = 0. 

As expected, the Poincare section has the same symmetry as the full trajectory, and the 
symmetry in each of the projections is the same: Since there are points on the attractor with 
trivial symmetry, T{A) = 1. The 10-fold rotational symmetry of the attractor indicates that 
S(A) = Z 5 X Z 2 = ZiQ. The reflections in D 5 x Z 2 = Dio act on the plane as reflections 
across lines at the angles 0, vr/10, 27r/10, ..., Ovr/lO. The attractor does not have any of 
these reflectional symmetries. The attractor in each Poincare section consists of 10 equally 
spaced points on a circle, and hence it has reflectional symmetry about 10 lines of reflection, 
but these are not at the “correct” angles and hence they are not symmetries of the attractor. 
(The limit cycles in the Vi and V 2 projections appear to have these “spurious” symmetries 
as well, perhaps because the (t)-waveforms are nearly sinusoidal.) 

By considering the time order that the points on the Poincare section are visited, we can 
determine that the attractor has the spatio-temporal symmetry Xj(t + T/10) = —Xj+i(t), 
where T is the period of the limit cycle. (Spatio-temporal symmetries are not included in 
S(A).) 
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Figure 5: A periodic attractor for 5 coupled van der Pol oscillators in the projections onto 
Vi, V 2 , Wi, and W 2 . The dots are where the trajectory intersects the invariant Poincare 
section wq = 0. Note that the symmetry is that same in each projection: only one is needed 
to observe the symmetry of the attractor. The Runge-Kutta step size is h = 0.01 throughout 
this paper. 
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Note that the projections of the limit cycle cross the lines of reflection symmetry (e.g. 
Im(vi) = 0) even though every point in the attractor has trivial isotropy. These are “non¬ 
generic” points that project to points with the wrong symmetry. This illustrates how the 
symmetry of the attractor can be too big. For example, consider the Poincare section in 
Vi- The 10 points of the attractor will shift around as a parameter is varied, and at isolated 
parameter values one of the points will just happen to fall on the line Im(vi) = 0, and 
we would conclude wrongly that = D 5 x Z 2 . Fortunately, such coincidences are very 

unlikely with chaotic attractors. 


3.3 Six coupled van der Pol oscillators 


As our second example we choose 6 oscillators because 6 is rich in factors. The smallest 
faithful representation of Dg x Z 2 is 3-dimensional: Vq © Vi, and the smallest representation 
isotropy equivalent to is 6 -dimensional: Vq © lA © F 2 © V3. 

First we explored Equation (|^) with 5 = 0. Although we did not do an exhaustive study, 
we never saw quasiperiodic or chaotic dynamics. Either the attractor was a single point at 
the origin, with symmetry T{A) = E(A) = Dg x Z 2 , or it was a limit cycle with nontrivial 
instantaneous symmetry, as shown in Table 1. The crucial hypothesis in Theorem |2.6| is that 


the instantaneous symmetry T{A) should be an isotropy subgroup of the observation space 
V. This is enough to proceed, however a short theoretical discussion describing the rest of 
Table 1 is illuminating. 

Assume that T{A) is nontrivial (that is, T{A) 7 ^ 1). To simplify notation, let T := T{A). 
Recall that A C FixM(T)- Hence S(FixM(T)) is the largest possible average symmetry that 
A can have. It is an interesting fact that S(FixM(T)) = A"r(T), the normalizer of T in F, 
which is dehned for the abstract group without reference to the representation. We dehne 
Tj’ := N-p{T)/T, which is the largest possible symmetry of A modulo the instantaneous 
symmetries. (Beware that 0 and @] use T' where we write F^.) Thus, the “effective” phase 
space is FixM(T), and the “effective” symmetry group is F^. 


T := T{A) 

V 

Ft 

FixR6(T) 

E(A) 

{p, k) = Dg 

Vo 

Z2 

{(a, a, a, a, a, a)} 

Dg X Z2 

{pa, k) = Dg 

V3 

Z2 

{(a, —a, a, —a, a, —a)} 

Dg X Z2 

= Z2 

V2 

D3 X Z2 

{(a, b, c, a, b, c)} 

(pw) — Zg X Z2 

(pV) = Z2 

v 

D3 X Z2 

{(a, —b, c, —a, b, —c)} 

(pw) — Zg X Z2 


Table 1: Observed symmetry of attractors of Equation (|3TD with n = 6 and 5 = 0. The 
hrst column is the instantaneous symmetry T, and V is the representation space for which 
T is an isotropy subgroup. The effective symmetry group and phase space are Tt and 
FixR 6 (T) © FixR 6 (T). The last column show the average symmetry E(A), which was always 
observed to be the same for a given T. 


Figure ^ shows the projection onto V 3 © IF 3 for a limit cycle with T{A) = {pa, k). The 
symmetry can be deduced by the projection onto V 3 alone, but since V 3 is one-dimensional 
we can draw pictures in V 3 ©hF 3 . Each oscillator is tt out of phase with its nearest neighbors. 
Within the effective phase space FixRi 2 (T) = R^, Equation (^I]) reduces to a single van der 
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Pol oscillator with effective symmetry P^ = Z 2 


X = y + 26y 

y = —X + y{{a — 2/3 + 27 — 25) — x^) 


(45) 


where ^3 = 6x and = 6y are the coordinates in V 3 ©hP 3 . The attractor has fnll symmetry, 
S(hl) = De X Z 2 , since the limit cycle has inversion symmetry. 



Fignre 6 : A periodic attractor of Eqnation (^I]) with n = 6, a = 1, (3 = —1, and 7 = 5 = 0. 
The symmetry is described in the second row of Table 1: T = {pa, k) . 

The cases where T = Z 2 rednce to Eqnation (|3lD with n = 3, so there are effectively 
3 van der Pol oscillators with D 3 x Z 2 symmetry. Fignre shows the Vi projection of a 
periodic attractor with T = {p^a), the last row in Table 1. We can determine T in two ways: 
First, the Xi coordinates of a single point on the trajectory match the pattern of FixR 6 (T) 
in Table 1. (The //* coordinates also match this pattern.) Alternatively, the signatnre of 
this instantaneous symmetry is that the projections onto Vq and V 2 are 0, while there are 
non-zero projections onto Vi and V3. Once the instantaneous symmetry is known, S can be 
determined by the Vi projection alone. We can think of F^ as {p^,K,a). The elements in 
Ft that give the symmetry observed in the Vi projection are (p^, a) = Z 3 x Z 2 = Zg. Thus, 
the full symmetry is E(A) = (p^, a, p^a) = (p, a). 

When 5 7 ^ 0 more exotic symmetries and quasiperiodic or chaotic dynamics are common. 
As an example, consider the one-parameter family with a = 1, /9 = 0, —2 > 7 > —5, and 
5 = —0.1. The instantaneous symmetry is T{A) = 1 throughout. The average symmetry 
and dynamics are summarized in Table 2. 

Unlike the cases with nontrivial instantaneous symmetry, we are forced to use a 3- 
dimensional projection, Vo©Ui, which is the smallest faithful representation space. We used 
the invariant Poincare section defined by tco = 0. We plotted two different 2-dimensional 
projections of Vq © Ui (the “top” view and the “side” view) on the computer screen. The 
top view is the Vi plane, and the side view shows the real part of vi horizontally and vq 
vertically. We used three colors to signify vq < —e, |no| < e, and vq > e, for an appropriate 
small e. Color figures of each type of attractor listed in Table 2 can be found at the web site 
[http : //odin.math. nau. edu/^jws/ . 

For this article, we show black and white figures of selected attractors. (The color would 
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Figure 7: A periodic attractor for 6 coupled van der Pol oscillators in the Vi projection. 
The dots are where the trajectory intersects the invariant Poincare section wq = 0. The 
parameters in Equation (^) are a = 1, /3 = 0.75, 7 = —0.5, and 5 = 0. The symmetry is 
described in the last row of Table 1: T = (p^cr). 


Dg XZ 2 

- 2.0 >Y > - 2.2 

Y =-2.3 Y =-2.4 

-2.5 >Y ^-4.2 -4.3 >Y^-4.5 

-4.6 > y >-5 

C 

Zg X Z 2 

Q 


C 


Ze 


Q 



Z 3 X Z 2 


P P 




Z 3 P 


Table 2: The observed dynamics and symmetry in Equation ( |^ for the family with 
n = 6, a = 1, /d = 0, 6 = —0.1 and 7 the parameter. The dynamics are either periodic 
(P), quasiperiodic (Q), or chaotic (C). The averaged symmetries are S(A) = Dg x Z 2 , 
Zg X Z 2 = (p, ct), Zg = (p), Z 3 X Z 2 = (p^,ct), or Z 3 = (p^). We changed 5 by steps of 0.1 
to produce the table, and found bistability only at 7 = —2.3. 
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be needed to distinguish between S(yl) = {p^,K,a) and S(A) = but we never 

observed either of these two symmetries.) 

The hgures we show in this article concern the bifurcation that occurs between 7 = —4.2 
and 7 = —4.3. From Table 2 we see that the symmetry changes from Z 3 to Zg x Z 2 . The 
number of conjugate attractors, which is |r|/|E| in general, goes from 8 to 2. Does this 
change of symmetry occur in 2 distinct bifurcations? Much to our surprise, the answer is no. 

As 7 decreases, the periodic orbit at 7 = —4.2 has a Hopf bifurcation to an invariant torus, 
and then the torus becomes crinkled, giving a presumably chaotic attractor at 7 = —4.2765, 
shown in Figure ^(a). From the Vi projection (the top view) there are two possibilities: 
E = (p^) or S = {p^,p^a). (Recall that p^a acts trivially on Vi.) The side projection in 
Figure ^a) shows that p^a ^ S. Hence the symmetry is S = (p^) = Z 3 , unchanged from 
the symmetry of the periodic orbit at 7 = —4.2 listed in Table 2. 




-20 Re(Vi) 20 -20 Re(v,) 20 

Figure 8 : Two projections of the periodic attractor when (a) 7 = —4.2765 and (b) 7 = 
—4.2766. The other parameters are a = 1, j3 = 0, and 6 = —0.1. The upper “top” view 
shows the projection onto Vi. In the lower “side” view, the horizontal coordinate is the real 
part of Vi and the vertical coordinate is Vq. Points on the invariant Poincare section are 
plotted whenever wq = 0. There are 50,000 points in(a) and 200,000 points in (b). 

Figure ^(b) shows the chaotic attractor at 7 = —4.2766. Here the average symmetry is 
E(A) = (p, (j) = Zg X Z 2 . It is clear from the side view that (p^,cr) is a subgroup of S(A). 
The attractor A in Figure ||(a), along with the conjugate attractors aA, p^A, and p^cA, 
merge into the single attractor in Figure |^(b) at a single bifurcation. There is no hysteresis, 
except for the fact that in going from (b) to (a) one of the four conjugate attractors is chosen 
randomly. The trajectory of the attractor in Figure H(b) spends several thousand iterates 


25 






near one the four “ghost” attractor left over from figure (a). Then it escapes and stays far 
from the ghosts for a few hundred iterates until it “falls” into any of the ghost attractors. 
Eventually, the trajectory visits all four ghosts. The symmetry of the attractor thus changes 
from (p^) to {p^,p^,a) = {p,cr) at a single bifurcation. 

This can be explained as a symmetry increasing crisis. Figure shows a blow-up of the 
boxed area in the side view of Figure ^(a). The folds in the attractor shown in Figure ^(a) 
indicate that the invariant circle has broken down in a standard bifurcation from quasiperi¬ 
odicity to chaos. We can infer the existence of a saddle-type periodic orbit of the Poincare 
map, indicated by an “x” in the inset of Figure ^(a). Somewhere between 7 = —4.2765 and 
7 = —4.2766 the attractor collides with the periodic orbit in a homoclinic bifurcation (the 
symmetry increasing crisis). At the bifurcation, one side of the one-dimensional unstable 
manifold is homoclinic, while the other side of the unstable manifold is in the basin of at¬ 
traction of all four ghost attractors. This is what causes the random hopping between the 
four ghost attractors. 




Figure 9: Blow-up of the two attractors in Figure H, showing the symmetry increasing crisis 
where 4 attractors merge into one. One million points on the invariant Poincare section are 
computed at each parameter value. The inset in (a) suggests that the dynamics is chaotic. 
A “ghost” attractor, and part of its conjugate under are seen in (b). 


This symmetry increasing crisis is different from the collision of two attractors discussed 
in 0 and [§. In a collision, the size of the attractor doubles. In the crisis we observed, 
the attractor more than quadruples in size. A similar bifurcation with two ghost attractors 
was described as an explosion in [j^. Crises of attractors, when no symmetry is present, 
are discussed in |]^. Typically, a crisis is accompanied by hysteresis, but in a symmetry 
increasing crisis there is no hysteresis. 


4 Conclusion 

We have given simple and easily verihable conditions to ensure that the projection of an 
attractor in an equivariant dynamical system has the same symmetry (generically) as the 
attractor. If the attractor has trivial isotropy, then the projection to a faithful representation 
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space is sufficient. This method is not always practical. For example, consider dynamical 
systems with permutation symmetry acting on R”. Then the smallest faithful represen¬ 
tation is (n — l)-dimensional. For a group like Sio acting on the projection technique is 
hopeless. On the other hand, symmetry detectives or any other method are also bound to 

fail. 

For group actions where the smallest faithful representation is 3-dimensional (or higher), 
progress can be made by looking at various 2-dimensional projections, using color, or by 
rotating a 3-dimensional projection of the object in real time on the computer screen. For 
example, with x Z 2 acting on R"^, the smallest faithful representation is three-dimensional. 
Figure ^ demonstrates how two different 2-dimensional projections, the top view and the side 
view, show the symmetry of the attractor. 
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